Photons from jet - plasma interaction in relativistic heavy ion collisions 
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We expose the role of collisional energy loss on high pr photon data measured by PHENIX collaboration 
by calculating photon yield in jet plasma interaction. The phase space distribution of the participating jet is 
dynamically evolved by solving Fokker-Planck equation. It is shown that the data is reasonably well reproduced 
when contributions from all the relevant sources are taken into account. Predictions at higher beam energies 
relevant for LHC experiment have been made. 



PACS numbers: 25.75.-q, 12.38.Mh 
Keywords: energy-loss, quark-gluon-plasma 



o 

(N 

< 

CN 



^ ■ 

Oh 1 
i ' 

43 



> 

CO 
CO 

o 

o 

oo 
o 



X 



I. INTRODUCTION 

Heavy ion collisions have received significant attention 
in recent years. Various possible probes have been stud- 
ied in order to detect the signature of QGP. Study of 
direct photon and dilcpton spectra emanating from hot 
and dense matter formed in ultra-relativistic heavy ion 
collisions is a field of considerable current interest. Elec- 
tromagnetic probes have been proposed to be one of the 
most promising tools to characterize the initial state of 
the collisions [l| . Because of the very nature of their in- 
teractions with the constituents of the system they tend 
to leave the system almost unscattered. In fact, pho- 
tons (dilepton as well) can be used to determine the ini- 
tial temperature, or equivalently the equilibration time. 
These are related to the final multiplicity of the pro- 
duced hadrons in relativistic heavy ion collisions (HIC). 
By comparing the initial temperature with the transition 
temperature from lattice QCD, one can infer whether 
Quark Gluon Plasma (QGP) is formed or not. 

There are various sources of photons from relativistic 
heavy ion collisions: (i) Direct photons are those which 
are produced in reactions of the type a b — > cj. One can 
subdivide this broad category of "direct photons" into 
"prompt", " prc-equilibrium" , "thermal" (from QGP as 
well as hadronic phase) and finally the "jet-photon" also 
called "jet conversion photon" (photons from jet-plasma 
interaction) depending on their origin, (ii) Decay pho- 
tons are basically the decay product of long lived sec- 
ondaries (ir° — > 77, 7] — > 77, p° — > 7T7T7, lo — > irj etc). 
The first calculation of thermal photons from quark mat- 
ter (QM) has been done in Ref. Q using hard thermal 
loop (HTL). It is shown that at a fixed temperature the 
contribution from QM is similar to that produced from 
hot hadronic matter. Here, our main concern is to calcu- 
late the jet-photons and compare them with the thermal 
photons from QM that has been calculated in Ref. @ . 

The jet conversion mechanism [|| occurs when a high 
energy jet interacts with the medium constituents via 
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annihilation and Compton processes. It might be noted 
that this phenomenon (for Compton process) has been 
illustrated quite some time ago [J] in the context of es- 
timating photons from equilibrating plasma, where, it is 
assumed that because of the larger cross-section, gluons 
equilibrate faster providing a heat bath to the incoming 
quark-jet. A comparison of the photons, calculated in the 
above scenario, (equivalent to photons from jet-plasma 
interaction) with the direct photons (thermal) shows that 
the former remains dominant for photons with pt upto 
6 GeV [1, HJ. It is to be noted that while evaluating 
jet-photon it is assumed in Ref. Q that the largest con- 
tribution to photons corresponds to p~ ( ~ Pqipq)- This 
implies that the annihilating quark (anti-quark) directly 
converts into a photon. Moreover, the quark jet might 
lose energy due to scattering with the constituents of the 
thermal bath before participating in Compton and anni- 
hilation processes. We include this effect on the photon 
productions. 

Collisional energy loss of heavy fermion has been cal- 
culated long ago in Ref. 0, Q using HTL approach. The 
same for the light quarks and gluons has been discussed 
in Ref. Q and revisited recently in Ref. However, its 
importance was demonstrated in the context of RHIC in 
Ref. Q . The measurements of non-photonic single elec- 
tron data [l^l show larger suppression than expected. 
These electrons mainly come from heavy quark decay 
where the radiative energy loss is suppressed due to dead 
cone effect. This observation has led to re-thinking the 
importance of collisional energy loss both for heavy as 
well as light quark [fl|. In view of this fact, collisional 
energy loss has been re-investigated in great detail (TlT - 
[T3 | in recent times. 

It is argued in Ref. [H that measurement of photons 
from such a novel process can provide direct informa- 
tion about the quark momentum. This is because of the 
assumption made in Ref. [|| that photons are predomi- 
nantly emitted at p 7 ~ p q . This implies that the thermal 
distribution of the participating parton is evaluated at 
the photon momentum. In this work we relax this as- 
sumption and calculate the photon yield from jet-plasma 
interaction. We consider photon production in the pt 
range 4 < pt < I4GeV. It is to be noted that to produce 
such photons the required energy of the participating jet 
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does not exceed (or remains below) the critical energy 
E c (-EWd = E co \\ at E = E c ) In fact, in this energy 
regime collisional loss seems to play important role [15l - 
Il7| . It is shown in Ref [H], [TtJ that the quenching factor 
for high pt hadrons can be accommodated within the 
framework of collisional energy loss only. 

Given the present scenario of energy loss mechanism 
in the context of RHIC data we, in this work, investigate 
the role of collisional energy loss as calculated in Refs. @, 
8:]. However, in order to see the effects of energy loss 
on jet-photon one should also incorporate the radiative 
energy loss and this has to be done in the same formalism 
in a realistic scenario. This has recently been done in 
Ref. [HI], where it has been shown that the neutral pion 
Pt spectra is sensitive to the inclusion of collisional and 
radiative energy loss. 

In the photon production rate (from jet-plasma inter- 
action) one of the collision partners is assumed to be 
in equilibrium and the other (the jet) is executing ran- 
dom motion in the heat bath provided by quarks (anti- 
quarks) and gluons. Furthermore, the interaction of the 
jet is dominated by small angle scattering. In such sce- 
nario the evolution of the jet phase space distribution is 
governed by Fokker- Planck (FP) equation where the col- 
lision integral is approximated by appropriately defined 
drag and diffusion coefficients. 

The plan of the paper is as follows. We give a brief 
description of photon production from QGP in section 
IIA. The evolution of jet quark and photon px distribu- 
tions are discussed in sections IIB and IIC respectively. 
We then briefly mention the necessary formulae for pho- 
ton production in initial hard collisions in section IID. 
Section III is devoted to the discussions of results and 
finally, we summarise in section IV. 



II. FORMALISM 
A. Thermal Photon Rate 

The lowest order processes for photon emission from 
QGP are the Compton scattering (q(q) g — > q(q) 7) and 
annihilation (qq — > gj) process. The total cross-section 
diverges in the limit t or u —> 0. These singularities have 
to be shielded by thermal effects in order to obtain in- 
frared safe calculations. It has been argued in Ref. [l9| 
that the intermediate quark acquires a thermal mass in 
the medium, whereas the hard thermal loop (HTL) ap- 
proach of Ref. [2(| shows that very soft modes are sup- 
pressed in a medium providing a natural cut-off k c ~ gT, 

We assume that the singularities can be shielded by 
the introduction of thermal masses for the participating 
partons. The differential cross-sections for Compton and 
annihilation processes are given by (2lj . 
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where m is the in-medium thermal quark mass, m 2 = 
"Zrrith 2 = 47ra s T 2 /3, a e and a s are the electromagnetic 
fine-structure constant and the strong interaction cou- 
pling constant, respectively. The static photon rate in 
1 + 2—^3 + 7 can be written as [l], Q 
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fi(Ei) , f 2 (E 2 ) and f 3 (E 3 ) are the distribution func- 
tions of 1st, 2nd and 3rd parton respectively, s, u and t 
are the usual Mandelstam variables. Mi represents the 
amplitude for Compton or annihilation process. Mi is 
the overall degeneracy factor. For Compton scattering 
Afi = 320/3 and for annihilation process Mi = 20 when 
summing over u and d quarks. 



B. Fokker - Planck Equation: Parton transverse 
momentum spectra 

As mentioned already in the introduction that the 
, 2 quark jet here is not in equilibrium. Therefore the cor- 
' responding distribution function that appears in Eq. J3]) 
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is calculated by solving the FP equation. The FP equa- 
tion, can be derived from Boltzmann kinetic equation 
(BKE)if one of the partner of the binary collisions is in 
thermal equilibrium and the collisions are dominated by 
the small angle scattering involving soft momentum ex- 
change [llB-lil. 

To arrive at the relevant FP equation from BKE we 
assume that there is no external force and therefore, 



d_ 
dt 



in+Vp-Vr /(p,x,f)=C[/( P> x ) t)] 



(4) 



Here, quarks have a phase space distribution which 
evolves in time and the collision term is evaluated by 
considering ultra-relativistic scattering of the quarks and 
gluons which eventually are expressed in terms of trans- 
port coefficients. For a Bjorken expansion [28| the prob- 
ability distribution is independent of the transverse co- 
ordinates and invariant under boosts in the z-direction 
and Eq. ^) is considerably simplified. Considering that 
the system expands in the longitudinal direction Eq. (U) 
takes the following form (2t| : 

a f (p,z,t) , v JJ^iA = c[f{ ^ t% (5) 



dt 



dz 



Here, v pz = p z /E p (for light partons E p = \p\). This 
equation can be simplified further for the central rapidity 
region which is boost invariant in rapidity, which implies 



f(PT,Pz,Z,t) = f(PT,P' z ,T). 



(6) 



Here p' z = y(p z — u z p), the transformation velocity u z = 
z/t, 7 = (1 — ul)^ 1 / 2 = t/r and r = \Jt 2 — z 2 denotes the 
proper time. Using the Lorentz transformation relation 
dr/dz\ z=0 = 0, 7 Z=0 = 1 and dp' z /dz\ z=0 = -p/t, one 
finds 



df Pz df 



dz 



t dp z 
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Therefore the Boltzmann equation takes the following 
form 



9f(PT,P z,t) 

dt 



\Pzt 



d_ _ Pz d 
dt t dp z 



f(pr,P„t) (8) 



Evidently in Eq. the second term on the left hand 
side represents the expansion while the right hand side 
characterizes the collisions. The latter can be written in 
terms of the differential collision rate Wp. q 

C\f(pr,P„t)] = J d 3 q[W p+q . tq f(p + q)-W p , q f(p)llO) 
which quantifies the rate of change of the quark momen- 



tum from p to p — q, W p 
represents scattering rates. 



<£T(p, q)/d 3 q, where T 



In a partonic plasma, small angle collisions, with para- 
metric dependence of 0(g 2 T), are more frequent than the 
large angle scattering rate. The latter goes as ~ 0(g 4 T). 
Therefore the distribution function does not change much 
over the mean time between two soft scatterings. This 
allows us to approximate /(p + q) — /(p)- In contrast, 
Wp,q; being sensitive to small momentum transfer, falls 
off very fast with increasing q. Therefore, we write 

^P +q . q /(p + q) - Wp,q/(p) + 3i^(W p , q /) 

+1 2^dBp-^ f) (u) 

With these approximation, Eq. [9] can be written as 
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arrive at the following equation 
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In Eq. ([TB"]) f(p,t) represents the non-equilibrium distri- 
bution of the partons under study, rj = (l/E)dE/dx, 
denotes drag coefficient, Bi = d((Ap z ) 2 ) /dt, B t = 
d((ApT) 2 ) /dt, represent diffusion constants along paral- 
lel and perpendicular directions of the propagating par- 
tons. 

The transport coefficients, rj, Bi and B t appeared in 
Eq. (|15p can be calculated from the following expressions: 
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(2tt) 5 J 2k2k'2p2p 
x (Af)^o/(|k|)[l±/(|k + q|)H 2 , 



-S(oj - v p • q)<S(w - v k • q) 
(17) 



in the small angle limit [3 |||. Here /(|k|,i) denotes the 
thermal distributions for the quarks (Fermi-Dirac) or glu- 
ons (Bose-Einstein) and stands for the statistical de- 
generacy factor for the i th parton species. The matrix 
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elements required to calculate the transport coefficients 
include diagrams involving exchange of massless gluons 
which render dE/dx and B^ t infrared divergent. Such 
divergences can naturally be cured by using HTL (20| 
corrected propagator for the gluons, i.e. the divergence 
is shielded by plasma effects. 

In the coulomb gauge the gluon propagator for the 
transverse and longitudinal modes are denoted by Dqq = 



A; and D; 



~q l q J /q 2 )A t with [3C 
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The HTL modified matrix element in the limit of small 
angle scattering takes the following form [8, 9] for all the 
partonic processes having dominant small angle contri- 
butions like qg — > qg, qq — > qq etc.: 

\M\ 2 = g 4 C R 16(E p E k ) 2 1 A, (©,,«) 

+ (v p x«).(v k xg)A f (? 0l9 )| 2 (20) 



Thus, to obtain px distribution of the jets at time t 
we need to convolve the Green's function with the initial 
spectrum [27j : 



(24) 



We use the initial partonjDT distributions (at the forma- 
tion time ti) taken from |a, l3l|: 



A'iVo 



d 2 p 0T dy ( 



1 2/o=0 
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POT \q ' 
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where AT is a phenomenological factor (~ 1.5 — 2) which 
takes into account the higher order effects. The values of 
the parameters are listed in Table. Q] 





RHIC 


LHC 


1 


q 


q 


q 


iV [l/GeV*\ 


5.0 x 10* 


1.3 x 10* 


1.4 x 10 4 


1.4 x 10 b 


P [GeV] 


1.6 


1.9 


0.61 


0.32 


a 


7.9 


8.9 


5.3 


5.2 



TABLE I: Parameters for initial parton pr distribution. 



where Cr is the appropriate color factor. With the 
screened interaction, the drag and diffusion constants can 
be calculated along the line of Ref. Q . For jet with en- 
ergy E » T Eqs.(5) and (6), in leading log approxima- 
tion, give (e.g. for q q — > qq) Q 



dE 
dx 

B, 
B, 



E 

IF? 



It 3 



(21) 



Similarly, drag and diffusion coefficients for the relevant 
processes can be calculated analogously. Having known 
the drag and diffusion, we solve the FP equation using 
Green's function techniques: If P(p, t\po,ti) is a solution 
to Ea. tfPo)) with the initial condition 



P(p,t = k\po,U) =6^(p-p ) 



(22) 



the full solution with an arbitrary initial condition can 
be obtained as 1241 



f(t,p)= / d 3 p p(p,t\p-b,u)f (p ) 



(23) 



where for the initial condition f(t = 0,p) = fo{po) and 
P(p, t\po, U) is the Green's function of the partial differ- 
ential Eq.(15). 



C. Space time evolution 

In this section we discuss how to obtain the space- 
time integrated rate of photons from jet plasma interac- 
tion using Bjorken hydrodynamics [28|. Note that, for 
jet photon, one of the distribution function appearing in 
Eq.(3) should be replaced by the phase-space distribution 
of the incoming jet. We assume invariant Bjorken corre- 
lation [32l | between the particle rapidity (y) and the space 
time rapidity (rj) to obtain the phase phase distribution 
of the jet. Now we have 



d 3 xd 3 p 
(2tt) 3 



f{x,p) =Ni 



(26) 



where Ni is the number of particles i and Vi is the spin- 
color degeneracy The phase-space distribution function 
for an incoming (quark) jet, assuming Bjorken n — y cor- 
relation [111], is as follows, 



fjet(r,p,U) 



(2ir) 3 V(r ± ) dN 

v q rpT d 2 prdy 
{2^) 3 r{r 1 _)Up T dN 
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where, ti is the formation time of the jet, and zq is its 
position in the QGP expansion direction. We consider 
the jets to be massless. It is also assumed that the jets 
do not change direction due to its interaction with the 



plasma particles. In such case fj e t can be factorizcd into 
a position space and a momentum space part and finally 
we obtain the phase space distribution of the jet at a 
later time if and at y = as(see Ref. [33| for details), 



fjet(r,p,t')\ v=0 



(2ir) 3 V{\w r \) U 1 



V q \Jt? - Zq 
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PT 
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( PT ,t')S(z ) (28) 



dN 
d' 2 pTdy 



(pT,t') can be obtained from Eg. ([24]) . v q is 



where 

the spin-color degeneracy factor for the incoming quark 
and 'P(luvl) is the initial jet production probability dis- 
tribution at the initial radial position w r in the plane 
zq = 0, where 

Kl = (r-(t'-U) Z-)-f 
\P\ 

= d \r cos 4>-t') 2 + r 2 sin 2 <j> for U ~ (29) 

and 4> is the angle in the plane zq = between the direc- 
tion of the photon and the position where this photon has 
been produced. We assume the plasma expands only lon- 
gitudinally. Thus using d A x — rdrdt' dtfidzQ we obtain the 
transverse momentum distribution (using Eqs.(|3]) and 
|)) of photon as follows: 



dN~< 
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= I d 4 x 
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16(2tt) 7 £ 7 
1 dN 
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dsdi\Mi\ 2 J dE 1 dE 2 

(PlT,t ) , 

VaE 2 2 + 2bE 2 + c 
(30) 



where ^prrdy (pit, t') can be calculated from Eq.p4)). 
f2i /3 are Fermi-Dirac or Bose-Einstein distributions. t c 
is the time when phase transition from quark matter 
(QM) to hadronic matter (HM) be gins and can be ob- 
tained by using Bjorkcn cooling law [28| . R is the trans- 
verse dimension of the system. <j) dependence occurs only 
in V{w r ). So the <fi integration can be done analytically 
as in Ref. [33l |. The temperature profile is taken from 
Ref. HI. 
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FIG. 1: Centrality fraction as a function of maximum impact 
parameter b m . 



so that perturbative expansion in powers of a s becomes 
meaningful. Thus the partonic cross-section, reaction 
rates etc. can be calculated without much ambiguity 
and with great degree of accuracy. Therefore, within this 
model the hard photons coming from initial hard parton- 
parton collisions can be estimated very accurately. 

In order to calculate reaction of the type h^hs — > 7 X 
(where Iia^b refer to hadrons), we assume that the en- 
ergy is such that the partonic degrees of freedom be- 
come relevant and they behave incoherently. The cross- 
section for this process can then be written in terms of 
elementary parton-parton cross-section multiplied by the 
partonic flux which depends on the parton distribution 
functions 34j . The energy scale (so-called factorization 
scale) for this to happen is denoted by Q 2 , the square 



of the momentum transfer of the reaction. Starting with 
two body scattering at the partonic level the differential 
cross-section for the reaction of above type can be written 
as M 



da 
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d 2 p T dy 



abc Jx * 
2 X a X b 



dx a G a/hA (x a , Q 2 ) G b/hB {x b , Q 2 ) 



da 



7T 2x a — xxe y dt 



(ab —} 7c). 



(31) 



where, xt = 2ptI\/s. The elementary partonic cross- 
sections for Compton scattering and annihilation process 
are given earlier. We also include photons from frag- 



D. Hard Photons 

The large p^-phenomenon in hadron-hadron collisions 
is well described by the perturbative QCD (pQCD) im- 
proved parton model. In this model it is assumed that the 
partonic structures of hadrons are revealed at high ener- 
gies. The strong coupling constant a s becomes 'weak' 



mentation process. This is accomplished by introducing 
the fragmentation function, Z? 7 / C (z, Q 2 ), when multiplied 
by dz gives the probability for obtaining a photon from 
parton c. Here z is fractional momentum carried by the 
photon. The differential cross-section is, therefore |35| . 

P^f = KJ2 f dx « f G a/hA (x a ,Q 2 ) 
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FIG. 2: 1Z as a function of maximum impact parameter b„ 
corresponding to RHIC and LHC. 



FIG. 4: Same as Fig.© at LHC. Tt = 0.897GeV and n 
0.073 fm/c. 
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where 
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2x b 

■.„ x a x T e~ 



2x a — XTe v 
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(33) 



2-x T e-y 

Once the photon production cross-section is obtained 




P T [GeV] 



FIG. 3: pr distribution of light flavors. The parameters (such 
as initial temperature T; and time n) correspond to RHIC 
experiment. T t = 0.446 GeV and n = 0.147 fm/c. 

from hadron-hadron collision we can now determine the 
direct photon production rates due to hard scattering 



between partons from nucleus-nucleus collisions at rela- 
tivistic energies. To do this we must note that the experi- 
mental data are given for a particular ccntrality. In order 
to take into account this experimental fact we introduce 
the centrality parameter (or the most inelastic fraction). 
It depends on the maximum impact parameter b m an can 
be calculated from the expression: 



f(bm) 



J^db(l-[l-T AB (b)a^ N ) AB ) 
J °°dh(l-[l-T AB (b)^ N ] AB ) 



(34) 



From this expression we extract b m relevant for a given 
experiment and use it to calculate photons from initial 
hard collisions and from parton fragmentation. Thus the 
yield becomes 



dN AB 
d 2 pr dy 

where 



(b m ) = TZ(b m ) 



da. 



7, hard 



da 



7, frag 



K{b m ) = (ABTab) = 



d?pT dy d 2 pT dy 



/ bm d 2 bABT AB (b) 



(35) 



f 6m db(l-[l-T AB (&)<7^]^) 



and 



T AB (b) - J d 2 sT A (s)T B (b 



(36) 



(37) 



is the nuclear overlap function. 



III. RESULTS 

First of all let us concentrate on the centrality mea- 
surements in an experiment. The photon measurement 
is done for a given centrality. For example, 10% cen- 
trality corresponds to f(b m ) ~ 0.1. In Fig. [T] we plot 
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FIG. 5: (Color online) pr distribution of photons at RHIC 
energies. The red (blue) curve denotes the photon yield from 
jet-plasma interaction with (without) energy loss. The solid 
(dotted) black curve corresponds to hard (thermal) photons. 
The magenta represents the total yield compared with the 
Phenix measurements of photon data [3^]. T; = 0.446 GeV 
and U = 0.147 fm/c 



the most inelastic fraction as a function of the maximum 
impact parameter b m for Pb-Pb collisions at LHC. We 
obtain b m ~ 5 fm for 10% most central collision for Pb- 
Pb system. Similar value is obtained for RHIC. 7Z vs b m 
given by Eq. I3U is plotted in Fig. [2] for Pb-Pb and Au- 
Au systems from which we obtain 7Z ~ 215(235) fm -2 at 
RHIC (LHC) for 10% centrality. We shall use these val- 
ues while estimating hard photon yields at RHIC (LHC) 
energies. We plot the transverse momentum distributions 
of quarks in Figs. [3] and 2] for different times at RHIC and 
LHC energies respectively where the initial distributions 
are taken from Eq. (|2"8j) . It is seen that as the time in- 
creases the quark stays longer in the medium losing more 
energy. As a result the depletion in the distribution func- 
tion is clearly revealed. It should be noted here that we 
do not include the induced gluon radiation which may 
further deplete the distribution at higher momenta. 

In Fig. [5] we show the pt distribution of photons from 
various processes which contribute at this high px range. 
It is observed that due to the inclusion of energy loss in 
the jet-plasma interaction the yield is depleted. Our cal- 
culation without energy loss is similar to that in Ref. Q 
at RHIC energies. However, at LHC energies, as we shall 
see below, the difference is by a factor of 2 — 3. Thus, we 
see that the assumption made in Ref could be valid 
for RHIC energies in the px range considered here, but 
at higher beam energies this is not a good approxima- 
tion. It is observed from Fig. [5] due to the inclusion of 
energy loss the rate is lowered by a factor ~ 1.5(1.7) 
at px = 4(14). This is more or less similar to what 
is obtained in Ref. [Hj|. The total photon yield con- 
sists of jet-photon, photons from initial hard collisions, 
jet-fragmentation and thermal photons. It is seen that 
Phenix photon data is well reproduced in our model. At 



high pt region the data is marginally reproduced. The 
reason behind this is the following. For high pt photon 
the incoming jet must have high energy where the radia- 
tive loss starts to dominate. Inclusion of this mechanism 
might lead to a better description of the data in high pt 
range. In order to understand the role of photons from 
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FIG. 6: (Color online) Same as Fig. ([5]) where in the total 
yield the contribution of photons from jet plasma interaction 
has been excluded. 



jet-plasma interaction in describing the Phenix high pt 
photon data we show in Fig. |6] where the contribution 
from jet-plasma interaction is excluded. It is seen that 
in order reproduce the data with 3 < pt < 6 GeV one 
must consider this extra source of photons. 
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FIG. 7: (Color online) pr distribution of photons at RHIC 
energies with T, = 0.350 GeV and n = 0.25 fm/c. The ma- 
genta line corresponds to total yield from all the sources as 
in Fig. ([5| with energy loss included in the jet-photon contri- 
bution. 
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To cover the uncertainties in the initial conditions for 
a given beam energy, we consider another set of initial 
conditions at a lower temperature T\ = 0.350 GeV and 
somewhat later initial time of Ti = 0.25 fm/c. The yield 
for this set is shown in Fig. [7] We see that the data is 
reproduced reasonably well. The inclusion of radiative 
energy loss will improve the situation further. 

As mentioned before, we also consider the high px pho- 
ton production at LHC energies. The contributions from 
various sources are shown in Fig. [7] Since the initial 
temperature in this case is higher, the plasma lives for 
longer time. Thus the energy loss suffered by the parton 
is more. As a result, the difference between the cases with 
and without energy loss is slightly more than what is ob- 
tained at RHIC. We do not consider strange quark as it 
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FIG. 8: (Colour online) Same as Fig. ([5]) at LHC energies. 
Ti = 0.897 GeV and n = 0.073 fm/c. 

is Boltzmann suppressed because of its high mass. How- 
ever, if one considers strange quark the effective number 
of flavours is taken as 2.5. For Np = 2(2.5) the lifetime 
of the plasma (t c —ti) is of the order of 1.9(2.6) fm/c (the 
values of T c are taken from Ref . [33] ) • So we see that the 
life of the plasma increases marginally. As a result mi- 
nor enhancement in the yield is expected. On the other 
hand, the degeneracy factor Mi is 20(30) Nj? = 2(2.5) in 
case of annihilation proceses. So the yield will be larger 
by a factor of the order of 1.5. Similar thing happens in 
the case of compton scattering. 



IV. SUMMARY 

We have calculated the transverse momentum distri- 
bution of photons from jet plasma interaction with colli- 
sional energy loss. It is shown that the assumption made 
in Ref. Q while calculating photons from jet - plasma 
interactions may not be good at LHC energies (the dif- 
ference is by a factor of 2 — 3). Phenix photon data 
have been contrasted with the present calculation and 
the data seem to have been reproduced well. We note 
that the yield of thermal photons as well as photons from 
jet-plasma interaction is very sensitive to initial temper- 
ature (Ti) time (ji) and the equation of state. It is shown 
that in order to describe the Phenix data in the domain 
4 < pt < 14 GeV the contribution from jet-plasma in- 
teraction is found to be important. The energy of the jet 
quark to produce photons in this range is such that colli- 
sional energy loss plays a dominant role here. In view of 
this fact we have considered collisional energy loss only. 
However, one should include both mechanisms in a con- 
sistent manner to describe the high px data beyond 14 
GeV or so. We note that the data is over-predicted un- 
less the energy loss is included. As we validate our model 
through the description of Phenix data we also predict 
the high px photon yield that might be expected in the 
future experiment at LHC. We observe that the contribu- 
tion from jet-plasma interaction is slightly more reduced 
as compared to the RHIC case. We notice that the in- 
clusion of the radiative energy loss will further reduce 
the yield at high px- That is expected to give a better 
description of the Phenix data in the high energy regime. 

Finally, it should be noted that the fragmentation pho- 
ton should, in principle, be calculated taking into account 
the energy loss of the fragmenting parton while traversing 
the plasma. This can be done by calculating the modi- 
fied pt distribution of the parton by using FP equation 
and then fragmenting it into photon as is done in case of 
high px hadron production in relativistic heavy ion colli- 
sions [TtJ • The yield will be depleted as compared to the 
result shown in the present work. This can lead to the 
situation where jet-photons and fragmentation photons 
may be comparable in the intermediate pt domain as is 
evident from the present results. However, in the present 
work we do not consider this aspect as our main concern 
is to see the effect of collisional energy loss on photons 
from jet plasma interaction. 



[1] J. Alam, S. Sarkar, P. Roy, T. Hatsuda and B. Sinha, 

Ann. Phys. 286, 159 (2000). 
[2] J. I. Kapusta, P. Litchard and D. Seibert, Phys. Rev. D 

44, 2774 (1991). 
[3] R. J. Fries, B. Muller, and D. K. Srivastava, Phys. Rev. 

Lett. 90, 132301 (2003). 
[4] P. Roy, J. Alam, S. Sarkar, B. Sinha, and S. Raha, Nucl. 

Phys.' A 624, 687 (1997). 



[5] E. Braaten and M. H. Thoma, Phys. Rev. D 44, 1298 
(1991). 

[6] E. Braaten and M. H. Thoma, Phys. Rev. D 44, 2625 
(1991). 

[7] M. H. Thoma, Phys. Lett. B273, 128 (1991). 

[8] A. K. Dutt-Mazumder, J. Alam, P. Roy, B. Sinha, Phys. 

Rev. D 71, 094016 (2005). 
[9] P. Roy, A. K. Dutt-Mazumder and J. Alam, Phys. Rev. 



C 73, 044911 (2006). 
[10] S. S. Adler et al., Phenix Collaboration, Phys. Rev. Lett. 

96, 032301 (2006). 
[11] A. Adil, M. Gyulassy, W. Horowitz and S. Wicks, Phy. 

Rev. C 75 044906 (2007); M. Djordjevic, Phys. Rev. C 74 

064907 (2006); T. Renk, Phys. Rev. C 76 064905 (2007). 
[12] S. Peigne, P. B. Gossiaux, and T. Gousset, J. High energy 

Phys. 04, 011 (2006). 
[13] J. Braun and H-J. Pirner, Phys. Rev. D 75, 054031 

(2007). 

[14] A. Ayala, J. Magnin, L. M. Montano, and E. Rojas, Phys. 

Rev. C77, 044904 (2008). 
[15] M. G. Mustafa and M. H. Thoma, Acta Phys. Hung. A 

22, 93 (2005). 
[16] A. Peshier, Phys. Rev. C 75, 034906 (2007). 
[17] P. Roy, J. Alam and A. K. Dutt-Mazumder, J. Phys. G 

35, 104047 (2008). 
[18] G. Y. Qin, J. Ruppert, C. Gale, S. Jeon, G. Moore, and 

M. G. Mustafa, Phys. Rev. Lett 100, 072301 (2008). 
[19] K. Kajantie and P. V. Russkanen Phys. Lett. B121, 352 

(1983). 

[20] R. D. Pisarski and E. Braaten, Nucl. Phys. B337, 569 
(1990); ibidNucl. Phys. B339, 310 (1990). 

[21] C. Y. Wong, Introduction to High Energy Heavy Ion Col- 
lisions, 1994, World Scientific, Singapore. 



[22] J. Alam, S. Raha and B. Sinha, Phys. Rev. Lett 73, 1895 
(1994). 

[23] B. Svetitsky, Phys. Rev. D 37, 2484 (1988). 
[24] G. D. Moore and D. Teaney, Phys. Rev. C 71, 064904 
(2005). 

[25] M. B. G. Ducati, V. P. Goncalves and L. F. Mackedanz, 

hep-ph/0506241. 
[26] J. Bjoraker and R. Venugopalan, Phys. Rev. C 63, 

024609 (2001). 

[27] H. V. Hees and R. Rapp, Phys. Rev. C 71 034907 (2005). 

[28] J. D. Bjorken, Phys. Rev. D 27 140 (1983). 

[29] G. Baym, Phys. Lett. B 138 18 (1984). 

[30] M. Le Bellac, Thermal Field Theory, Cambridge Univer- 
sity Press, Cambridge, 1996. 

[31] B. Muller, Phys. Rev. C 67 061901R (2003). 

[32] Z. Lin and M. Gyulassy, Phys. Rev. C 51, 2177 (1995). 

[33] S. Turbide, C. Gale, S. Jeon and G. D. Moore, Phys. Rev. 
C 72 014906 (2005). 

[34] J. Pumplin, D. R. Stump, J.Huston, H. L. Lai, P. Nadol- 
sky, W. K. Tung, J. High Energy Phys. 012 0207 (2002). 

[35] J. F. Owens, Rev. Mod. Phys. 59 465 (1987). 

[36] S. S. Adler et al., Phys. Rev. Lett. 98 012002 (2007) . 

[37] J. Braun and H. Pirner, Phys. Rev. D 75 054031 (2007). 



